Temporal dynamics in the one-dimensional quantum Zakharov equations for plasmas 
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Abstract 

The temporal dynamics of the quantum Zakharov equations (QZEs) in one spatial dimension, which describes the nonlinear 
interaction of quantum Langmuir waves (QLWs) and quantum ion-acoustic waves (QIAWs) is revisited by considering their 
solution as a superposition of three interacting wave modes in Fourier space. Previous results in the literature are modified 
and rectified. Periodic, chaotic as well as hyperchaotic behaviors of the Fourier-mode amplitudes are identified by the analysis 
of Lyapunov exponent spectra and the power spectrum. The periodic route to chaos is explained through an one-parameter 
bifurcation analysis. The system is shown to be destabilized via a supercritical Hopf-bifurcation. The adiabatic limits of the 
fully spatio-temporal and reduced systems are compared from the viewpoint of integrability properties. 

PACS numbers: 52.25.Gj; 05.45.Mt; 52.35.Mw. 
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I. INTRODUCTION 



The Zakharov equations (ZEs) are one of the most 
important models in plasma physics community [l|, in 
which high-frequency Langmuir waves (LWs) are nonlin- 
early coupled with the low-frequency ion-acoustic waves 
(IAWs). In this context, the quantum Zakharov equa- 
tions (QZEs) are the modified version of the classical 
ZEs, including a quantum correction associated with the 
Bohm potential Such QZEs are also deduced from a 
multiple time-scale technique applied to a set of quantum 
hydrodynamic (QHD) equations under the quasineutral- 
ity assumption. Recent works (see e.g. Refs. [sHf!) in- 
dicate that much attention has been paid to investigate 
the dynamics of such one-dimensional (ID) QZEs. Very 
recently, a more comprehensive work on the dynamics of 
LWs has been studied by Haas and Shukla in a three- 
dimensional quantum Zakharov system Q- The arrest 
of Langmuir wave collapse by quantum effects, predicted 
by a variational approach in Ref . [7j , was later confirmed 
with rigorous estimates and systematic asymptotic ex- 
pansions by Simpson et al @. 

In contrast to classical ZEs, the QZEs can not be re- 
duced, in the adiabatic limit, to a nonlinear Schrddinger 
equation (NLSE) Q. Rather, it follows a coupled system 
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for the envelope electric field and the density fluctuation, 
whose complete integrability is not assured. However, 
the system can be decoupled in the adiabatic as well 
as semiclassical limit 0, Q, where the solitons can be 
found to be more stable than in fully degenerate cases 
0. Other recent developments on the QZEs can be 
found in the literature 0-0, • Marklund Q studied 
the kinetic theory of LWs interacting with quantum ion- 
acoustic waves (QIAWs), where it was shown that the 
combined effects of partial coherence and quantum cor- 
rection tend to enhance the modulational instability (MI) 
growth rate [3j. The temporal dynamics of the QZEs has 
been studied by Haas by means of a variational approach 
[J] . It has been found that the quantum coupling param- 
eter plays a destabilizing role on localized structures, or 
Langmuir wave packets. More mathematical treatments 
on the QZEs are the analysis of the underlying Lie sym- 
metry group and the derivation of some exact solu- 
tions |1(J, LLL| ■ A Galerkin-type approximation was used 
by Misra et al. [f| to reduce the QZEs to a set of ordinary 
differential equations (ODEs) for the temporal dynam- 
ics. This system was shown to exhibit hyperchaos (more 
than one positive Lyapunov exponent). However, while 
qualitatively correct, unfortunately the reduced model 
considered in this work was flawed by some algebraic in- 
consistencies. 

The primary goal of the present work is to revisit 
the temporal behavior of the QZEs as a superposition 
of three interacting wave modes in Fourier space and to 
rectify the previous results [5[- In addition, the evidence 
for the existence of periodic limit cycles, chaotic as well 
as hyperchaotic attractors of the Fourier-mode ampli- 
tudes are presented through the analysis of bifurcation 
diagram, power spectra as well as the Lyapunov expo- 
nents. There are relevant differences of the respective 
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periodic, chaotic and hyperchaotic regimes in parame- 
ter space. Our findings thus extend both qualitatively 
and quantitatively the previous results as well as exhibit 
some new features not reported in the earlier investiga- 
tions [3, [al- 



ii. SIMPLIFIED MODEL 



The ID QZEs read @ 



.dE d 2 E 



dt dx 2 
d 2 n d 2 n 
dt 2 dx 2 



H- 



dx A 
dx 1 



= nE, 

d 2 \E\ 2 
dx 2 



(1) 
(2) 



where E = E(x,t) is the envelope electric field and 
n = n(x,t) is the plasma density fluctuation (measured 
from its equilibrium value) . In Eqs. (JTHUl , the same set of 
dimensionless quantities of Ref. [2| is employed. In par- 
ticular, H = hui/ kb T e is a quantum coupling parameter 
expressing the ratio between the ion plasmon energy and 
the electron thermal energy, where H is Planck's constant 
divided by 2-7T, is the ion plasma frequency, kb is the 
Boltzmann constant and T e is the electron fluid tempera- 
ture. The formal classical limit H — > yields the original 
Zakharov system 

The finite-dimensional temporal dynamics of the QZEs 
was studied in Ref. In this investigation, the ap- 

pearance of hyperchaos was established by means of the 
derivation of two positive Lyapunov exponents and the 
analysis of the Kaplan- Yorke dimension. The dynamics 
of this (low-dimensional) simplified system was conside- 
red for a wide parameter range of the system, including 
the quantum coupling parameter H. However, in the 
derivation of this simplified model the change of variables 
in Eqs. (7) and (8) of this work is not fully consistent. 
Indeed, Eq. (3b) of this paper implies d = —Ani sin</>, 
while Eq. (3c) would imply d = — 2ni sin</>, using the 
notations of Ref. Q. The reason for the contradiction 
is that the transformations [Eqs. (7), (8) of Q] do not 
respect the conservation of plasmon number as given by 
Eq. (5) of the cited paper. Similar difficulties also ap- 
pear in the literature about the classical Zakharov system 
[I2I HH, in the derivation of reduced ordinary differen- 
tial equations simulating the classical (H — 0) temporal 
as well as spatio-temporal system. Hence, in this Sec- 
tion we will rederive the basic set of equations assuming 
a truncated expansion scheme and produce the correct 
simplified model. The resulting system is then shown nu- 
merically to be compatible with the existence of chaotic 
as well as hyperchaotic attractors, in qualitative agree- 
ment with the conclusions in Ref. [f§]. In order to access 
the low-dimensional dynamical behavior derivable from 
Eqs. (p~H5]), we follow Refs. [5| and [13] considering a 
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FIG. 1: Hamiltonian function H in Eq. fT9J), showing distinct 
qualitative properties for different values of the parameter /x: 
H = "v/OlS (upper panel), /i = a/3 (lower panel). 



few-modes expansion through the ansatz 

E = y/N sin exp(z6» ) + V2N cos (|j cos(fca;) exp(i6>i) 

(3) 

n = N + rti cos(kx) , (4) 

where N is the conserved plasmon number, k is the fixed 
wave- number of the excitation and a — a(t), 6q — 9o(t), 
9\ = 6\ (t) and n\ = n\ (t) are real time-dependent quan- 
tities. It is important to notice that Eq. (|3]) implies 
\E\ 2 = 7V+oscillatory terms, in agreement with the con- 
servation of plasmon number. 

Inserting Eqs. (J3HH) into Eqs. ((T}]21) and computing 
the zero wave-number term and the coefficient of the 
term proportional to cos(kx) yields four equations, cor- 
responding to the real and imaginary parts. Further, 
the higher-order harmonics proportional to cos(2/cx) are 
disregarded, in the spirit of the present few modes ap- 
proximation. The real and the imaginary parts of the 
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independent term gives 



where 



= -V2m 



sm tp , 



V2 



7 = - N - m cos <^3 cot (^j 



(5) 
(G) 



where the 'dot' represents derivative with respect to time 
t and 



tp = Vo 



(7) 



Now, considering the real part of the coefficient of the 
term proportional to cos(fcir) we get 

01 = -N - k 2 (l + H 2 k 2 ) - ^nicos^ tan . (8) 

The imaginary part of the coefficient or the term propor- 
tional to cos(fca;) provides no further information. The 
combination of Eqs. ([6] [8j) gives 



= k 2 (l + H 2 k 2 ) - \pln\ castp cot a. 



(9) 



To close the system, it is necessary to substitute the 
ansatz of Eqs. (J5H11) also into Eq. ©. Ignoring again 
higher-order harmonics yields 

hi + k 2 (l + H 2 k 2 ) ni = -V2Nk 2 costp sina. (10) 

The system of Eqs. ([S]),© and ([TO")) provides a finite- 
dimensional closed set of ordinary differential equations 
for a, tp and n\, whose properties will be analyzed nu- 
merically in the next Section. 

Of particular interest is the case where the density fluc- 
tuations respond adiabatically to excitations. In this si- 
tuation the second derivative term in Eq. (fT0|) can be 
neglected, yielding 



rii = 



y/2N 

1 + H 2 k 2 



cos tp sm a . 



(11) 



Inserting this last result into Eqs. ([5]) and ©, we get the 
two-dimensional dynamical system 



N 



1 + H 2 k 2 



sm2ip sina , 



tp = k 2 (l + H 2 k 2 ) + 



2N 



1 + H 2 k 2 



(12) 

cos 2 tp cos a . (13) 



Actually, Eqs. (|12D and (|T3")) contain only one relevant 
free parameter fx, as can be proved by introducing the 
rescaled time variable 



Nt 



1 + H 2 k 2 



so that 



da/dr — sin 2^ sina, 
dtp/dr = /i 2 + 2 cos 2 tp cos a , 



(14) 



(15) 
(16) 



k(l + H 2 k 2 ) 



(17) 



In this formulation, the role of quantum effects, contained 
in the free parameter H, are hidden in the rescaled time 
r and also in fi. Thus, the quantum effects tend to slow 
down the dynamics, since a larger H causes dr/dt to be- 
come smaller in magnitude. Moreover, it follows from Eq. 
((TBI that the fixed points can exist only when // < y2- 
But, since \x is a monotonically increasing function of H, 
the conclusion is that quantum effects tend to suppress 
the existence of equilibrium points. 

In addition, Eqs. (|15H16[) form a completely integrable 
dynamical system, as can be better seen in writing it in 
the generalized Hamiltonian form 



da 



J 



dU 

dtp ' 



dtp 



-J 



dU 

da 



(18) 



where J = esc a and the first integral H playing the role 
of Hamiltonian function is 



% = jJL cos a 



cos 2 tp sin 2 a . 



(19) 



The system (jT8")) can be written in canonical Hamiltonian 
form using Darboux coordinates (u,v) = (— cosa, tp), so 
that du/dr = dH/dv,dv/dr = —dH/du. Most impor- 
tantly, Eq. (fTS)) is manifestly non-chaotic. Notice, how- 
ever, that the level surfaces of the constant of motion H 
are not compact, and thereby violating a necessary con- 
dition for Liouville integrability. In addition, the Poisson 
structure is singular in the sense that J is not well-defined 
for a — lir, where I is an integer. Figure 1 exhibits some 
typical graphs of H for /i = \/Q,5 and [i = \J% showing 
distinct qualitative properties of the constant of motion. 

On the other hand, the adiabatic limits of the simpli- 
fied finite-dimensional system and the original infinite- 
dimensional system are to be compared. For the infinite- 
dimensional system, the adiabatic limit is get disregard- 
ing the second-order time derivative of the density fluc- 
tuation in Eq. The resulting equations are 



.dE d 2 E 



dt dx 2 



\E\ 2 E = H 2 



d*E 



dx 2 



n=\E\ 



E W)> (20) 
(21) 



which in the formal classical limit (H — > 0) reduce to 
the usual NLSE, which is well-known to be completely 
integrable. However, in the quantum case, the adiabatic 
limit still shows a coupled nonlinear system, whose prop- 
erties are not yet completely known. In contrast, the 
simplified dynamics associated to the present few-modes 
ansatz has shown to be integrable. Notice that the sys- 
tem (l2"0H2"l"j) can be decoupled taking both the semiclas- 
sical (H -C 1) as well as adiabatic limits, so that the 
substitution n = — \E\ 2 is allowed in the right-hand side 
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of Eq. (j2"0")) . In this case one obtains 

^ + a^ + |s|£; - F U^ -i5 -^J- (22) 

Equation (|22l) can be used as the starting point for study- 
ing quantum perturbations of the classical NLS soliton 
solutions jg|. 

III. NUMERICAL RESULTS 

The system of Eqs. ((5J), ([9j and (fT0|) can be recast as 

±1 = — V2 x 3 sinx2 , (23) 

x 2 = k 2 (1 4 H 2 k 2 ) - V2x 3 cosx 2 cotxi , (24) 

x 3 = X4 , (25) 

id = -k 2 (1 4 H 2 k 2 ) x 3 - V2Nk 2 cosx 2 sin x x , (26) 

where, for convenience, we have redefined the variables 
as a = x\, (p = X2, fix — x 3 , h\ = x±. 

Eqs.([23Hini) were numerically solved by using the sixth 
order Runge-Kutta-Fehlberg scheme with the step length 
h = 0.01 and initial values as x\ = 0.1, X2 = 0.2, 
Xs = 0.3, a?4 = 0.1. A Monte-Carlo search on the param- 
eter space was conducted to find the possible dynamics 
of the system, namely periodic, chaotic or hyperchaotic 
regimes. In order to establish the existence of irregular 
dynamics we have investigated the Lyapunov exponents 
spectra [l4[ . Chaos or hyperchaos is characterized by the 
presence of one or two positive largest Lyapunov expo- 
nents respectively. It is sufficient to calculate only the 
three largest Lyapunov exponents. These largest Lya- 
punov exponents were calculated by integrating Eqs. (|23I 
|26| in order to have average estimates of them over the 
attractors [l5j]. Figure 2 shows examples of two largest 



Lyapunov exponents with respect to the quantum pa- 
rameter H and for constant N = 1.5 and k = 0.8 such 
that /i < \[2 is satisfied. We observe that the system is 
chaotic (one positive exponent) for a small range of H, 
while it is hyperchaotic (two positive exponents) for a 
wide range of values of H . In the classical case, i.e., for 
H = 0, we also find two positive Lyapunov exponents 
indicating the hyperchaotic features even in the classi- 
cal ZEs [l| as well. It is to be noted that the system 
may experience some other features at the points where 
some peaks of the first (maximum) Lyapunov exponent 
seem to be correlated with the depletions of the second 
exponent. However, since the variations of the exponents 
are shown with respect to a system parameter, namely 
H, there may be a possibility for the existence of chaos 
rather than the hyperchaotic orbits at those points, es- 
pecially where the second Lyapunov has negative peaks. 
Figure 3 shows an example of the hyperchaotic dynamics 
of the system in the x\ — x 3 plane for H = 0.2, N = 1.5 
and k = 0.8. It seems to show an intermittent switching 
between two regions of the phase space, one to the left 
and other to the right of x\ ~ 3.2. These can be veri- 
fied from the corresponding time series for x\ and x 3 as 
shown in Fig. 4. Here the upper panel clearly explains 
this intermittency near x\ = 3.2, and both the time series 
show the aperiodic nature of the system which are very 
common in the context of chaotic dynamics. The chaotic 
dynamics of the system can be established by calculat- 
ing the power spectrum corresponding to the variable x 3 . 
Here one can measure the power spectrum as the square 
of the modulus of the complex Fourier coefficients corre- 
sponding to x 3 with frequency as the inverse of the period 
of the signal [lg|. However, measuring the amplitude of 
the pulse gives rise a slight modification of the power of 
the same obtained by using the Fast Fourier Transform 
(FFT). The latter is used to compute the discrete Fourier 
Transform (DFT) of the signal, as well as the magnitude 
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FIG. 4: The time series corresponding to the variables x\ 
(upper panel) and X3 (lower panel). The parameter values 
are the same as in Fig. 3. 

and phase of the transformed signal [13] • The 'abs' func- 
tion is used to obtain the magnitude of the data and the 
'angle' function to obtain the phase information [18j . and 
unwrap in order to remove the phase jumps greater than 
7r to their 2tt complement. These are illustrated in Fig. 
5. This figure clearly shows the chaotic dynamics of the 
waveforms with low- frequency and high amplitudes, as 
well as exhibits the broadband spectral features of the 
system. The route to chaos from a single periodic orbit 
can be explained by means of a one-parameter bifurca- 
tion analysis in the domains of dynamical variables. We 
have considered one fixed value of H = 0.2, two different 
values of k, namely k — 0.9, 0.95 and allowed the plasmon 
number N to vary so as to satisfy the inequality /i < \/2. 



FIG. 5: The amplitude (absolute value, see the upper panel) 
and the phase (lower panel) of the hyperchaotic signal cor- 
responding to the variable x$. The parameter values are the 
same as in Fig. 3. 



The three-dimensional pictorial views of the bifurcation 
diagram with x±, X3 and N are shown in Fig. 6 for dis- 
tinct values of k. These explain that the system loses its 
stability through a supercritical Hopf-bifurcation (hb), 
which gives the birth to periodic limit cycles at N « 0.78 
for H = 0.2, k = 0.9 (upper panel) and N m 0.6 for 
H = 0.2, k = 0.95 (lower panel). The nature of the 
supercritical Hopf-bifurcation can be seen in the X1-X2 
plane from Fig. 7 for N — 1.5, H = 0.2 and different 
k: (i) k — 0.3 (upper panel) and (ii) k = 0.5 (lower 
panel). Notice that in the case of Hopf-bifurcation, we 
have generated a program with varying k from k = 0.1 to 
1.0. Figure 7 shows illustrations of two different values 
of k. The trajectories in both the upper and lower panel 
of Fig. 7 are generated with different initial conditions, 
and finally they converge to attracting or repelling limit 
cycles. In both the cases the initial values are taken from 
the unstable sets. For values of k larger than k = 0.5 no 
more attracting or repelling limit cycle is found to coex- 
ist. From Fig. 7 one can also observe that a limit cycle 
(the blue curves or e.g., the orbits appeared first from the 
bottom at X\ < and at x\ > ) is surrounded by the 
unstable equilibriums (red colored curves). To further 
examine the chaotic features, one can also calculate the 
Fourier transform of the signal [l(| corresponding to the 
variable x\. Figure 8 represents the distribution of the 
Fourier coefficients in the complex plane corresponding 
to the chaotic signal x\iov H = 0.2, N — 1.5 and k = 0.8. 
It shows how wide and dense are the Fourier coefficients 
in the complex plane. This distribution could be useful 
to quantify chaos from numerical points of view. How- 
ever, the detail analysis is beyond the scope of the present 
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periodic orbit is analyzed by means of one-parameter 
bifurcation analysis. Moreover, the reduced temporal 
dynamics is shown to be integrable in the adiabatic 
limit. This is in contrast to the adiabatic limit of the 
spatiotemporal dynamics described by the QZEs whose 
integrable properties deserves further study; whereas 
the finite-dimensional reduced system is shown to evolve 
into periodic, chaotic as well as hyperchaotic orbits 
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FIG. 6: The snap shots of an animated movie showing the 
three-dimensional views of the bifurcation diagrams with vari- 
ation of the parameter N and for a fixed H = 0.2. Two dif- 
ferent values of k are: k = 0.9 (upper panel) and k = 0.95 
(lower panel). Different colors may correspond to the pixels 
of different projections on the "floor" . 
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IV. CONCLUSION 



FIG. 7: Supercritical Hopf- bifurcation in the x\ —X2 plane for 
N — 1.5, H = 0.2 and different k: (i) k = 0.3 (upper panel) 
(ii) k = 0.5 (lower panel). 



The nonlinear interaction of quantum Langmuir waves 
and quantum ion-acoustic waves is analyzed in terms 
of a superposition of three interacting wave modes 
in Fourier space. Previous works on both classical 
[l2l Il3| and quantum Q Zakharov equations have been 
rectified and modified. The chaotic behaviors of the 
reduced temporal system have been identified by the 
analysis of Lyapunov exponent spectra as well as by the 
analysis of power spectrum. The hyperchaos has been 
characterized by the presence of two positive Lyapunov 
exponents. Also, the route to chaos from a single 



To conclude, we believe that the results presented in 
this work would be helpful for better understanding the 
salient features of such ID QZEs. 
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